Quantum measurement backaction from a BEC coupled to a mechanical oscillator 
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We study theoretically the dynamics of a a hybrid optomechanical system consisting of a macro- 
scopic mechanical membrane magnetically coupled to a spinor Bose-Einstein condensate via a nano- 
magnet attached at the membrane center. We demonstrate that this coupling permits us to monitor 
indirectly the center-of-mass position of the membrane via measurements of the spin of the con- 
densed atoms. These measurements normally induce a significant backaction on the membrane 
motion, which we quantify for the cases of thermal and coherent initial states of the membrane. We 
discuss the possibility of measuring that quantum backaction via repeated measurements. We also 
investigate the potential to generate non-classical states of the membrane, in particular Schrodinger 
cat states, via such repeated measurements. 

PACS numbers: 42.50.Lc, 42.50.Wk, 42.79.Gn, 07.10.Cm 



I. INTRODUCTION 

While the foundations of Quantum Measurement The- 
ory were already laid down in the early days of quantum 
mechanics, quantitative studies are a relatively recent de- 
velopment [H [2] ■ These studies are now driving signifi- 
cant advances in quantum information science and quan- 
tum metrology, with major progress initiated by the de- 
sire to test Bell's inequalities on the one hand, and by the 
quest for gravitational wave detection on the other (3J |4j . 
A common aspect of these and related studies is the need 
to quantify, control, and possibly exploit the quantum 
backaction of one or a series of measurements on a quan- 
tum mechanical system. 

An important development in this context is the 
rapid progress witnessed by cavity optomechanics, which 
makes it increasingly realistic to consider the use of 
mechanical systems operating in the quantum regime 
to make precise and accurate measurements of feeble 
forces and fields pT. In many cases, these measurements 
amount to the detection of exceedingly small displace- 
ments. In this context, hybrid systems consisting of cou- 
pled atomic (or molecular) and nanomechanical systems 
may prove particularly useful. The robust and scalable 
infrastructure provided by NEMS/MEMS devices, cou- 
pled with the high precision measurement capability of 
quantum gases [l2lU4j . makes them an attractive com- 
bination for sensitive force measurements, as well as for 
a quantitative study of dissipation and decoherence pro- 
cesses at the quantum-classical interface. As a result, 
there are ongoing experimental [5j |6] and theoretical j7j- 
[TTj efforts toward coupling mechanical systems to atomic 
ensembles. 

The system that we consider in this paper consists of a 
mechanical membrane magnetically coupled to a spinor 
Bose-Einstein condensate, an arrangement somewhat 
similar to a system previously considered in Refs. [7lH0], 
The magnetic coupling, generated via a nano-magnet an- 




FIG. 1: Proposed experimental setup for backaction detec- 
tion, involving a magnetic particle located at the center of 
a vibrating membrane and polarized along the z-axis. This 
setup produces a spatially inhomogeneous magnetic field that 
is detected by the BEC. The BEC's long axis is along the 
z-direction, the probe light (red arrow) is along y, and the 
membrane oscillations are along the a;-axis. 



chored on the membrane, entangles the membrane to the 
spin state of the BEC. We may then make an indirect ob- 
servation of the center-of-mass position of the membrane 
through a measurement of the spins of the condensed 
atoms. While such magnetic coupling and concomitant 
entanglement is possible between the membrane and a 
thermal gas of spinor atoms (or even a single atom) , the 
BEC has two advantages over the alternatives when used 
as a position sensor. Namely, bosonic stimulation be- 
tween the atoms leads to increased signal strenth, and 
two-body collisions are suppressed, giving rise to a much 
longer phase coherence time when compared to thermal- 
ized atoms. Now, the indirect position measurement can 
be achieved generally speaking in two ways: through a 
strong, or projective, measurement of the spin state of 
the BEC, or through a weak, or dispersive, measurement 
of the spins. An example of the former would be through 
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a Stern-Gerlach-type measurement. Here, the atoms are 
subjected to a spatially inhomogenous field seperating 
out the different spin components. The different spin 
populations can then be measured via a standard absorp- 
tive imaging technique. A likely candidate for the latter 
type of measurement would be observation of the Larmor 
frequency through a phase contrast imaging technique, 
as described in Ref. [12] . Since the Larmor frequency is 
proportional to the local magnetic field - which is mod- 
ulated by the motion of the magnetic domain attached 
to the oscillating membrane - it provides a measure of 
the membrane motion. Either type of measurement of 
the BEC spin can induce a backaction on the membrane, 
modifying its position and/or momentum in proportion 
to the strength of the measurement. The main goal of 
this project is thus to quantify the effect of using the BEC 
as a position sensor on the membrane and to evaluate the 
feasibility of measuring the backaction of such a measure- 
ment. In the case where the backaction is small we can 
indeed use the BEC to make ultraprecise measurements 
of the membrane, and when it is large, it may provide a 
means to observe the effects of a quantum measurement 
on a macroscopic object. So far such a quantum effect 
has only been observed in dilute and isolated systems 
such as ultracold gases [16] . 

Due to the scope of this endeavor, we find it convenient 
to divide our results into three papers. In this, the first 
paper, we present the experimental setup and develop 
a simplified Hamiltonian to describe the interaction of 
the membrane and condensate. We derive a general for- 
mulation of the system state after a single or multiple 
arbitrary measurements of the condensate's spin state, 
and then present results for specific cases where the mea- 
surements are projective. We show that in principle the 
membrane can be put into a non-classical state by such a 
procedure. In the second paper we extend these consid- 
erations to weak measurements, focusing on quantifying 
the backaction of indirect observation of the Larmor fre- 
quency via phase contrast imaging. In the third paper we 
will exploit the backaction of repeated measurements and 
coherent control to achieve quantum control of the state 
of the membrane, allowing in principle the production of 
squeezed or Schrodinger cat states. Because the effects 
of environmental coupling, i.e. decoherence and dissipa- 
tion, are critical to coherent control protocols, we will 
also include these effects in the computations of paper 3. 

This paper is organized as follows; Section II intro- 
duces the system under consideration and derives a model 
Hamiltonian that describes the magnetic coupling of the 
membrane to the condensate. Section III discusses the in- 
direct measurement of the membrane center-of-mass po- 
sition via a strong measurement of the spin of the conden- 
sate, assuming that the membrane is initially in a thermal 
state. In particular, we present the post-measurement 
Wigner function of the membrane and interpret its main 
features in terms of a backaction parameter. Section IV 
then turns to the case of repeated measurements on this 
system, illustrating how a sequence of measurements pro- 



vides a direct signature of the backaction. Section V dis- 
cusses the possibility of producing highly non-classical 
states of the membrane, considering specifically the case 
where it is initially in a coherent state. We show that 
repeated measurements typically leave the membrane 
in a non-classical state characterized by a non-positive 
Wigner distribution. Finally, Section VI is a summary 
and outlook. Some technical details of the calculations 
and experiment are included as appendices. 

II. MODEL 

As mentioned in the introduction, the system under 
consideration consists of a micromcchanical membrane 
of fundamental mode frequency ui m and effective mass 
m, whose center-of-mass is oscillating, perhaps under the 
influence of a weak force. Our goal is to develop and an- 
alyze a quantum measurement scheme that permits us 
to characterize the center-of-mass motion of that mem- 
brane. 

The scheme that we envisage involves integrating this 
membrane into a hybrid system whose other component 
is an elongated spinor Bose-Einstein condensate with 
long axis z, see Fig. 1. The condensate is subjected to 
a static magnetic field B$ along the quantization axis z. 
The membrane is magnetically coupled to the conden- 
sate via a magnetic domain anchored at the center of the 
membrane. 

We assume for simplicity that the magnetic domain is 
a point dipole located at the origin and polarized along 
the z-axis, /i m = /i m z . The magnetic field due to that 
dipole is 

B(rH^4( 3 ^m-?)?-Mm). (!) 

The small spatial variation of the resulting magnetic field 
along the long axis of the condensate results in a variation 
of its Larmor precession frequency. As we show in the 
following, this dependence permits us to characterize the 
expectation value of center-of-mass mode of oscillation of 
the membrane. 

The Zeeman interaction between a single atom in the 
condensate and the total magnetic field Bo + B(r) is de- 
scribed by the interaction Hamiltonian 

V = - M a-(Bo+B(r)) 

= fisgp (F X B X + F y B y + F Z (B + B z )) , (2) 

where Ft is the i fh component of the spin-1 operator, [is 
is the Bohr magneton, and gp = 2 for alkali atoms. Be- 
cause of the extremely low temperature of the condensate 
and the relatively small trapping frequencies employed 
(order 10 - 100 Hz), we will neglect the kinetic energy of 
the atoms and the trapping potential in our model. As 
a result of the spatial dependence of the magnetic field, 
atoms at different positions along the long axis of the con- 
densate precess at different frequencies and dephase over 
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a period of time. It is this z-dependent phase difference 
that is picked up in spin population measurements. 

We note that the transverse dependence of the mag- 
netic field also results in an inhomogeneous broadening of 
the Larmor frequency. The effects of this broadening can 
be minimized, the more elongated along z and tightly 
confined in the transverse directions the condensate is. 
In the following we consider for simplicity a condensate 
that is almost one-dimensional and confined to a region 



close to y — 0, so that B y 



and F y B y 



0. Further- 



more, close to z — 0, i.e. for the fraction of the BEC 
directly above the dipole, the magnetic field B(r) is pre- 
dominantly in the z- direction provided that xq, the equi- 
librium distance between the BEC and the membrane, is 
much greater than the relevant coordinates y and z. In 
that case we can ignore the effects of B x and B y alto- 
gether. As shown in detail in Appendix A, for small 
displacements x m of the membrane compared to xq, the 
single atom magnetic coupling Hamiltonian ([2| reduces 
then to 



is the membrane Hamiltonian, and we have dropped (and 
will continue to drop) the subscript in membrane dis- 
placement for compactness, x rn — > x. By explicitly com- 
bining both terms of H and completing the square for x, 
we express the total system Hamlitonian as 
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H = + -m^ n (x + AF Z ) + {HNn L0 F z - hSflF*} 

(8) 

where 



A = fiBgpNB'Jmoj 2 



(9) 



which we call the backaction parameter in anticipation of 
the following sections, 



^lo = HBgFB c /h, 



(10) 



which is simply the Larmor precession frequency of the 
atomic spins, and 



SQ = mLo z m A 2 /2h 



(11) 



V = fi B gFF z 



Bn 



4ttx% 



(-x + 3z m ) 



The first and second terms of this expression are inde- 
pendent of time, while the third term, proportional to 
x m , varies sinusoidally in time. We exploit this property 
by rewriting the magnetic Hamiltonian as 



V = vb9fF z (B c + B' v x ri 



where 



B c 
BL 



= Bn- 



jl jl rn 

AttxI 1 



Airx* 



(5) 



By treating the condensate as N non-interacting spin- 
1 atoms under the mean-field approximation and by as- 
suming that the magnetic field acting on the atoms in the 
detector region is roughly constant, we can easily arrive 
at the many-atom magnetic interaction Hamiltonian. It 
is simply 

= N^ B g F F z {B c + B' v x m ), (6) 

where *b(x) is the annihilation operator for an atom at 
position x (i.e. it is a Schrodinger field operator). The 
total system Hamiltonian of the hybrid BEC-membrane 
system is then 



Here 



H — H rn + Vbec- 



H m = p /2m + mw m x /2 



(3) III. SPIN MEASUREMENT 

In this and the following sections, we evaluate the back- 
action on the membrane of an arbitrary measurement op- 
erator that acts on the spin degree of freedom of the BEC. 
We consider the cases of a single and multiple spin mea- 
surements for a membrane initially prepared either in a 

(4) thermal or in a coherent state. The analytical expressions 
provided for the post measurement density matrix (and 
subsequently the Wigner function) are for a general mea- 
surement operator. However, to illustrate the backaction 
effects in an intuitive manner, we consider in this paper 
the particular case of simple projective measurements of 
F y , corresponding to all atoms being in the F y = 0, +1 
or —I state. While this is not the most likely outcome of 
a typical Stern-Gerlach measurement, it does serve as a 
simple demonstration of backaction induced dynamics. 

The measurement proceeds in the following way: at 
time t = a tt/2 pulse is applied to the condensate, 
preparing all atoms in \F X = 1) state. Following that 
preparation stage, the atomic spins precess about the z- 
axis until a later time t\ when the spin population mea- 
surement is made. Additional measurements can be per- 
formed at later time intervals t i: i — 2,3,.... Further 
precise details of the experimental setup are presented in 
Appendix D. 

Equation Q already gives a clear indication of the 
backaction of the measuring apparatus - the condensate - 
on the membrane. As a result of their coupling the mem- 
brane Hamiltonian is modified from being a harmonic os- 
cillator centered at the origin to one that is shifted by the 
quantity AF Z , indicating that the backaction depends on 
(7) the outcome of a specific spin measurement. Here, we 
give an explicit description of the measurement process 
by evaluating the pre- and post-measurement density op- 
erator of the membrane, and the corresponding Wigner 
function. 
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A. Density operator 

We assume that the membrane and the condensate are 
initially uncorrelated, 

P = Pm{0) ® /Obec(0) 

and denote the initial density matrix elements of the BEC 
as p a p, where p a ^ = (ct\pBEc\P) , and a, /3 are the vari- 



ous spin states, a,/3 = {0, ±1}. We also assume for now 
that the membrane center-of-mass is initially in thermal 
equilibrium with temperature T. Defining 

Tj = huj m /2k B T, 

its density matrix elements in position space are then 



p m (xf,Xi,t = 0) = (xf\p m \xi) 



mu>„ 



tanh r\ x exp 



mu)„ 



((xf + Xi) 2 tanh 7/ + (xf — Xi) 2 coth 77) 



(12) 



For t > the spin components of the condensate undergo a Larmor precession about the z-axis. Since [Hbec, H v 
V] = 0, we can use the Baker-Hausdorff relation to re-express the propagator U(t) = exp(—iHt/h) as 



U(t) 



_ e -itH B Ec/h e -it(H m +V)/h 



This allows us to find the evolution of the system density matrix in a straightforward way. After an interaction time 
ti, this evolution results in the matrix elements of the density operator of the membrane + condensate system to 
become 



,Xf\p(t)\/3,Xi) = p° a/3 exp[-iQ L a(a - 0)t x + iSQ L (a 2 - £ 2 )*i]y 



muj m 



tanh?7 x 



(13) 



exp 



rnuJm 

4h 



I (a;/ + Xi + (a + fi)A{l — cosw m ii)) 2 tanh 77 + (xf — X{ + (a — /3)A(1 — cosuj m ti)) 2 coth 77 



+4iA sin u m ti(ctXf — (3xi) + 2iA 2 (a 2 — j3 2 ) sinw m ^i(2 — coso; m ii)}] . 

I 



That is, the interaction of the membrane with the BEC 
displaces its center-of-mass motion in both position and 
momentum by amounts that depend explicitly on the 
spin components a and f3, as well as on the backaction 
parameter A. 



B. Single measurement 



As already discussed a measurement of arbitrary type, 
but dependent on F y , is carried out on the BEC at time 
t\. The post-measurement density matrix of the mem- 
brane depends on the measurement outcome, and is given 
by 



(Xf\Pn 



J_tr BEC (WV ys ). (14) 



Here <p is the outcome of the F y measurement, P(4>) is 
the probability of that outcome, is the Kraus op- 
erator corresponding to the effects of the measurement 
on the BEC's quantum state, and p sys is the complete 
system density matrix given in (13). Though Eq. (14) 



is true for arbitrary measurements (i.e. for arbitrary 
Kraus representations), in later plots and numerics, we 
take the operators to be projectors onto the eigen- 
states of F y , namely, the 3 operators \F y = j)(F y =7], 
7 = {0, ±1}. This simplified situation displays the salient 
features of measurement back-action without requiring 
extensive computations. Such a scenario could be real- 
ized through e.g. a Stern-Gerlach aparatus. As men- 
tioned above, we will return to the precise computation 
of operators corresponding to the planned experi- 
mental setup (Appendix D) in subsequent work. 



C. Phase Space Representation 

The effect of the measurement on the state of the mem- 
brane can be visualized particularly clearly in terms of 
its Wigner distribution function 



W(x,p) = —J dte-Wix + £/2\p\x - e/2) (15) 
For a harmonic oscillator in a thermal state, we have 
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W(x, p, t = 0) = — tanh 77 exp 
irn 



mto m / 2 



tanh 77 



an expression that should be contrasted to the post-measurement Wigner function, which is found to be 



W(x,p,h) = ^tanhr?-^y^W| Q p^exp[-in i0 (a- P)h+i5n L (a 2 - p 2 )h] 



(16) 



(17) 



exp 



muj r . 



x+-(a + /3)(l-cos(u; m t 1 ))) 

2 ) \ mu>. 



+ —(a + (3) sin(w m /ji) ) [tanh// 



+iA(a — j3) xsin(w m ti) 



— — (l - cos(cj m ii) ]+^(a + P) sin(w m <i) 
mui m \ 12 



V 



At time t\ the first line of Eq. (17) simply describes the 
imposition, due to the evolution of the BEC, of a phase 
on each term of W(x,p) which is dependent on the mea- 
surement outcome and spin indices. More interesting are 
the last two lines in that expression: The second line de- 
scribes a term-by-term shift in the initial Gaussian prob- 
ability distribution The new phase-space center of each 
term in the Wigner function depends on the spin indices 
of that term, and it rotates through phase space at the 
membrane frequency. 

As the third line in W(x,p,t{) is imaginary, it results 
in interference between the terms of the Wigner function 
(visible below in Fig. [2]). From inspection of Eq. (17 1, 
we see that the interference will increase with larger A. 
The physical significance of this interference can be un- 
derstood from the fact that the initial state of the BEC, 
\F X = 1), is not an eigenstate of the interaction Hamilto- 
nian (proportional to F z ). So, one can think of the con- 
densate as experiencing three interaction Hamiltonians 
simultaneously, one for each of its spin components, and 
it is the interference between them that leads to the os- 
cillations. We remark that the oscillations in the Wigner 
function can be seen as long as 



.4 



- Z j) 



y tanh ? 



(18) 



where x 



zp 



action parameter of Ao = NA sa = 9 x 10 m (for com- 
parison, A — 0.22x zp ). Fig. [2 shows the resulting post- 
measurement Wigner functionfor this specific value of Aq 
(figure c), and for backaction parameters of 0.1A (figure 
b), and O.OlAo (figure a), assuming that a measurement 
of F y with result 1 was made after an interaction time of 
<i = n/uj m . 



As alluded to by Eq. (18), in order to observe the os- 



cillations in the post-measurement Wigner function, we 
can either increase the temperature or increase A. How- 
ever, increasing the temperature leads to dissipation and 
decoherence losses that are ignored in the present anal- 
ysis, but result of course in a fast thermalization and 
associated smoothing of W(x,p). As mentioned above, 
these effects will be discussed in future work. A more 
promising approach to observe quantum interference ef- 
fects is to increase A (Eq. (|9|), either by increasing the 
number of atoms in the effective detection zone, or by 
increasing B' v via a decrease of xo. Since B' v scales as 
1/xq, this may be the easiest way to reach the regime 
of observable Wigner function oscillations. Note however 
that for decreasing xq the simplified interaction Hamilto- 
nian ([6| becomes less accurate as the components of the 
magnetic field along x and y become more important, so 
the approximations made in Eq. ^ will not be as valid. 

The expectation value of the center-of-mass position 
(x) of a membrane in thermal equilibrium is zero, and its 
variance is 



D. Backaction 

To illustrate the effect of backaction, we assume the 
following membrane parameters for the remainder of the 
paper (unless stated otherwise): uj m = 27rl0 6 rad/s, 
m = 5 x 1CT 13 kg, /j m = 2x 10~ n A • m 2 , and an initial 
temperature of 4K. The static external magnetic field is 
B a =0.1 Gauss, and the condensate is x a — 5 x 10~ 6 
m away from the membrane, resulting in a single-atom 
backaction parameter ^4 sa = 9 x 10~ 21 m. We assume 
that N = 10 5 atoms experience the same magnetic field 
in the detection region, yielding then an effective back- 



a 2 (x) = (x 2 ) - (x) 2 = (x 2 ) = — — cothry. (19) 

2muj m 

Immediately following the measurement, the membrane 
is no longer in thermal equilibrium, and (x) in gen- 
eral. For large backaction parameters the oscillations 
in W(x,p) become quite significant. In that regime the 
BEC is a poor position sensor, since its coupling to the 
membrane significantly perturbs the outcome of subse- 
quent measurements, see Fig [2] Here the measurement 
creates a significant change in the phase-space distribu- 
tion of the membrane that invalidates any 
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FIG. 2: (Color online) Wigner distribution function of the 
membrane after one measurement giving the result F v = 1, for 
several values of the backaction parameter: (a) A = O.OlAo: 
(b) A = O.lAo and (c) A = A , with A = 0.22^ 'h/2muj m - 
The evolution time is t — ir/w m in all three cases. The Wigner 
function has been multiplied by 10 6 to make the axes legible. 



information gained about the position. However, as we 
show in Section V, a high value of A is beneficial for 
preparing the membrane in highly non-classical states 
via repeated measurements. As with all schemes for 
state preparation involving repeated measurements, its 
not very efficient for highly excited initial states, as is 
the case in Fig. [2]- At 4 K, and for the membrane pa- 
rameters of this example, the mean phonon occupation 
number is 8.3 x 10 4 . 

To investigate the efficiency of our setup as a position 
sensor, it is therefore appropriate to consider the limit of 
small Aq. Consider for concreteness the specific example 
where the outcome of the spin measurement is F y = 1. 
Ignoring then terms of order A 2 , and for t\ — it/£Ilq, we 
find 



IV. SUCCESSIVE MEASUREMENTS 

Generally, successive measurements of the conden- 
sate's spin result in an accumulation of backaction ef- 
fects in the membrane. At the same time, they offer the 
potential for the coherent control of the center-of-mass 
motion of the membrane. To address such situations, we 
now consider the effect of a succession of measurements 
on the state of the membrane. The similarity between 



Eqs. (12 1 and (13 1 suggests that it should be possible to 



find a closed form for the density matrix of the compos- 
ite system after an arbitrary number of measurements 
on the BEC. Indeed, this is the case and such a form is 
presented below. 



(x) = -— sm((jj m TT/Q L0 ) coth?7, (20) 

with (a; 2 ) remaining constant to lowest order in A. We 
then have 

h A 2 

o- 2 (x) pm = coth?7 — sin 2 (w TO 7r/n L0 ) coth 2 77 

2mLO m 4 

(21) 

where the subscript "pm" indicates post-measurement. 
The minimum backaction occurs for 77 — > 00, or T — > 0. It 
also vanishes when the membrane frequency is an integer 
multiple of Larmor precession frequency, 

uj m = nfl L0 , (22) 

where n in an integer, in which case it is possible to carry 
out stroboscopic QND measurements of the membrane 
position. 



A. Analytical results 

We first consider the situation where the membrane be- 
gins in a thermal state, and assume that the interaction 
between the membrane and the BEC, (initially prepared 
in the F x — 1 state, is turned on at time t = 0. After 
time tx, a first measurement is performed on the spin of 
the BEC. The probe-BEC interaction time t mea s required 
to carry out that measurement is taken to be negligible 
compared to the other characteristic times of the system 
- in practice this is true for phase contrast imaging. The 
direct effect of the measurement is only on the BEC, and 
can be formally described by a Krauss operator that 
depends explicitly on the outcome of the measurement. 
Following that first measurement the system evolves uni- 
tarily for an additional time t2, becoming re-entangled. 
A second measurement is then performed, acting on the 
BEC with an operator , and this process is repeated 
n times. The set of operators are related to the above 
W<t> by = M^M*. 

Immediately following the nth measurement, the ele- 



ments of the BEC-membrane system's density matrix are given by (see Appendix B for more details) 
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(a,x f \p(t)\P,Xi) 



ai,... a n ;fii,...f) n 



x exp 



x exp 



n 



E ki^miai - ft) - - ft 2 )} 



-i 

i=i 
■muj, 



tanh?7 



["^T ^ Xf + X ^ A + Xl + A[ftA) 2 tanh?? + (x f + X[a]A - x» - X[(3]A) 2 cothr? 



+ 4iA(P[a](s/ + X[a]A) - P[/3]{x t + X[/3]A)) + 2iA 2 {(j)[a] - </>[ft)} 



(23) 



where 

M {i) 
X[a] 



(<t\M®\t), 

n 

Cn - E (°"' _ 0"i-l)cOsT i!n , 



i = l 



P[cr] = 



E - o"i_i) sin T i)ri , 
»=i 

E( fJ i - 0-j-i) 2 sinTj, n cos T i)n 



(24) 
(25) 

(26) 



/ = 1 

n— 1 n 



- 2 E E ( (7 ~ 'i-l)( Cr J Sin T J> C0S r i,n)( 27 ) 

i—1 



E tk if « < j and = otherwise. (28) 



The sums on OLi and ft run over {— 1,0, +1}, and in 
Eqs. (25) - (27) we have ao = ft = wherever they 



appear. To be clear, in the equations above, the quan- 
tities ti are the intervals between measurements rather 
than the times themselves. Also, care should be taken as 
A, P, and <f> are functions of the spin indices and thus 
must be recomputed for each term in the sum. Note also 
that the density matrix (23) is unnormalized. Its trace 



is equal to the probability of the particular sequence of 
measurement outcomes described by the specific set of 



operators {AfW} occurring in it. 

A few remarks are in order before turning to a discus- 
sion of numerical results: First, the density matrix of the 
membrane can be obtained from Eq. ( 23 ) by a partial 



trace. It is the sum of many different contributions from 
various shifted (in position and momentum) thermal en- 
sembles. A (numerical) computation problem does arise, 
though, as calculating the density matrix for n measure- 
ments requires summing over 9 n terms. In our numer- 
ics we are thus restricted to few-measurement scenarios. 
The timing of the measurements is also very important 
to the BEC functioning as a detector. In one extreme 
case, if the measurements are made exactly at the nat- 
ural frequency of the membrane, the interaction will be 
completely masked, except for the tiny second-order ef- 
fect of 5D,l- This supports the claim that position is very 
nearly a stroboscopic QND variable. 



The structures of Eqs. (23) and (25)-(27) indicate that 



the results of the earliest measurements continue to be 
as important as those of later measurements. This is 
because in the absence of dissipation, there is no atten- 
uation of the information gained nor of the backaction 
induced by the measurements. Including the effects of 
thermal dissipation reduces and eventually erases this 
memory effect. 



The Wigner function (15) of the membrane after 



measurements can be derived from ( 23 ) and is found to 
be 



W(r n) - 1 V M (n) M( n_1 ) M (1) o « M (1)t M (n ~ 1)f M (n)t 

x exp ^-i^ti(n£ (Q!i-A)-^i(Q!i ^ tanh? ? 



x exp 



x + ^(X[a] + X\p])) 2 + f-f- + £(P[a] + P[/5 

/V "Mum & 



tanh 77 



+ iA(P[a]-P[/3]) ( 2: 

2 / ' /"-',„ 



^(0H-^])| , (29) 



r 



a form that clearly illustrates the shift in position and 
momentum of the various terms in the sum. 



B. Numerical results 

One simple way to detect the effect of quantum back- 
action on the state of the membrane is to consider a se- 
quence of 2 measurements carried out in succession at 
times t\ and £2, and to compare the outcome of the last 
measurement to the outcome of a measurement at that 
same time but skipping the first measurement at t\. 
As a concrete example, we consider, for the parameters 
of the previous section, the following two scenarios: 

1. F y is measured once, at t<x — Tt/u> m (i.e., half the 
oscillator period); 



2. F y is measured first at t± = 7r/2a; m and then at 

< 2 = Tt/u m . 



By simply calculating the trace of Eq. (23) for the 
3 possible operators M^'s in scenario 1 we obtain the 
probabilities 



P(F y = ±1) 

P{Fy = 0) 



0.375 
0.25 



(30) 



In the second scenario 9 traces need to be evaluated in- 
stead of just 3. Summing then the three terms corre- 
sponding to the same final value of F y , so as to obtain 
the total probability of obtaining a particular value for 
the second measurement, we find 



P(F y = ±1) = 0.344 
P(F y = 0) = 0.312 



(31) 



The simple example shows that backaction of the inter- 
mediate measurement should be readily observable, since 
it changes the probability for the three possible outcomes 
of the measurements at t 2 by a significant amount, of the 
order of several percent. 

A slightly more complete look into the effects of inter- 
mediate measurements is provided by the Wigner func- 
tions shown in Figures [2] and [3j In Fig. [2] we show a series 
of plots for the case of a single measurement at t = 7r/w m , 



for different values of A (O.OlAo, O.lAo, and Aq, respec- 
tively). In Fig. [3] we increase the number of measure- 
ments from 1 to 2. The time of the final measurement is 
the same as in Fig. [2j 7r/o; m , but an intermediate mea- 
surement is made at time t\ = 7r/2w m . The outcome of 
the final spin measurement is again F y — 1, but, for the 
sake of direct comparison with Fig. |2j we average over all 
possible results of the intermediate measurement. This 
is equivalent to an experiment in which the outcomes of 
the intermediate measurement are discarded or ignored. 
As expected, the Wigner functions where an intermedi- 
ate measurement are noticeably different than the one 
lacking an intermediate measurement, hinting at a sim- 
ple way to characterize the impact of quantum backaction 
on the membrane dynamics. 



V. INITIAL COHERENT STATE 

We now turn to the situation where the center-of-mass 
state of the membrane is a coherent state, an initial con- 
dition that can be prepared by driving the membrane 
with a classical force. The main result of this section is 
to demonstrate that after n measurements, the center-of- 
mass of membrane is split into a superposition of up to 
3" discrete coherent states, hinting at the possibility of 
generating a macrsocopic Schrodinger cat. 



A. Evolution of the Coherent State 

The assumption that the membrane is initially in a 
coherent state allows us to eschew the density matrix 
formalism for the moment, and the initial state of the 
composite system is given by 

|*(0)) = D(a + ib )\0 

mem >®$> Q |a) (32) 

cc 

where J2a \ c a\ 2 = 1 anc ^ 



D(u) = cxp(ua^ — u* a) 



(33) 
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is the displacement operator for the center-of-mass state nal state for the hybrid membrane-BEC system (see Ap- 
of the membrane. Following a procedure similar to that pendix C) 
for the initially thermal membrane, we arrive at the fi- 




FIG. 3: (Color online) Wigner distribution function of the 
membrane in the case of two measurements: the final spin 
measurement, at time n/uj m , yields the outcome F v = 1 . 
As in Fig. j2|, the three Wigner functions are plotted for 
increasing values of the backaction parameter, A, and the final 
Wigner functions are averaged over all possible outcomes of 
the intermediate measurement. 



/ 

|*(0)) = ^M("^...M« QiCctl exp (n L0 a % ~5n L af 



ck,cki ,...,a n 



with 

a[a] = clq cosTi in + bo sinTi^ 
b[a] = b cos Tx >n - a sin Ti >n 



er 



A 
A 



X[a], 



= J2 ai 



A 



a (sinT M - sinT M _i) - fo (cosT M - cosT M _i) + V](<7j - crj-i)(smT jti - sinT i:i _i) 



i AQ\a\ \ + lb[a])lQmem) g |a) (34) 



(35) 
(36) 

(37) 



j'=i 



and all other definitions are as before. This wave-function 
is not normalized; the probability of obtaining a par- 
ticular sequence of measurement outcomes is given by 

(*(*)!*(*)>• 

Because the real part of the displacement is propor- 
tional to (x) and the imaginary part is proportional to 
(p), we have a very clear picture of the physics arising 
in this situation. The membrane's initial displacement 
oscillates back and forth semi-classically. However, each 
time the BEC's spin is measured, the membrane's wave 
function splits into three distinct components, each hav- 
ing received a different kick from its interaction with 
the different possible spin orientations of the BEC. This 
would indicate that the membrane was indeed put into 



a cat state, except for the fact that the kick is propor- 
tional to A, which is typically at most comparable in size 
to the zero point oscillations x zp . If, going forward, we 
can increase A in an experimentally realizable setting, 
this system may provide an excellent demonstration of a 
macroscopic object put into a highly non-classical state. 



B. Wigner function 

Following a sequence of measurements, the Wigner 
function of the initially coherent state of the membrane 
motion becomes 
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W(x,p) 



1 \ " M (n) M (n-1) A/ r(l) A/fWt M^ 1 ^ 

exp ( -z^t i (O i0 (a l - A) ~ 5n L {aj - 



x exp 



mu>„ 



, (x - x zp (a[a] + a\j3] + i{b[a] - b[(i]))f + ( - x zp {b[a] + b[(i] - i{a[a] - a[/3])) 

ft \ ///_•„, 



.4 



+ i(a\0\b[a] - a[a]b[/3}) - i~ — (Q[a] - Q[ 



l -(a[a]-a[P\Y - l -{b[a]-b[f3]f 



(38) 



Figure [4] shows the post-measurement Wigner function 
after successive measurements. The key point here is that 
repeated measurements can lead to very non-classical 
states, as is evident from the resulting negative valued 
Wigner functions. As can be seen in Eqs. (34) and (38 1, 



each measurement introduces different phase factors to 
the initial coherent state, along with splitting it into dif- 
ferent coherent states. It is the quantum interference be- 
tween the different coherent states thus generated that 
gives rise to the non-classical Wigner functions. Rela- 
tively few measurements are required to generate such 
non-classical states. In fact, it is possible to turn a co- 
herent state into a non-classical state after only one mea- 
surement. Fig. [4] shows one such case. Here, ao = bo = 1 
and the result of the first measurement after a time in- 
terval of t\ = Tr/uJ m is F y = 0. Repeated measurements 
after equal time intervals lead to states that resemble 
displaced Fock states and other more complex states. 

If the initial displacement is small in magnitude (as 
in Fig. [4|, production of this non-classical behavior is 
rather insensitive to the measurement timing. For this 
particular case, with the outcome F y = 0, it is seen for 
t\ in the approximate range 0.67r/w m — 1.87r/w m . On 
the other hand if the initial displacement is large (for 
example, ao = 50\/2, &o = 0), the non-positive Wigner 
function is much harder to find; in this case it only oc- 
curs from tx m 0.77w/u m — 0.79%/ui m . Nevertheless, in 
either case, successive measurements can build up some 
very interesting non-classical states. It may be possible 
to obtain a high degree of probabilistic quantum state 
control of the membrane if the dissipative effects do not 
wash out the interference too quickly. 



VI. CONCLUSION 

In conclusion, we have demonstrated that by coupling 
a magnetic membrane to a spinor BEC, we can monitor 
and manipulate the position of the membrane. A pro- 
jective measurement procedure induces significant back- 
action that can be measured for reasonable experimental 
parameters. We investigated the effect of this interaction 



for different initial membrane states, namely thermal and 
coherent states. We discussed the possibility to measure 
backaction of a quantum measurement on the membrane 
via repeated measurements and the potential to generate 
cat states of the oscillator via such repeated measure- 
ments in the case of an initial coherent state. This would 
be a major accomplishment as such a state has not yet 
been prepared for a solid, macroscopic object. 

As mentioned above, we will in future work look at 
the effects of using a more dispersive (and more exper- 
imentally plausible) measurement scheme based on ob- 
servation of the Larmor precession. Additionally, we will 
study the effects of including dissipation on the system's 
dynamics. Other possibilities for extending the model 
even further include inclusion of a general coupling to 
other spin components, i.e. an interaction Hamiltonian 
of the form F • B, because at short distances the x and y 
gradients of the magnetic field become significant. Also, 
the BEC is a spatially finite system, so we may wish to 
exploit the ability to measure multiple "pixels" of the 
condensate in order to gain better information about the 
membrane's position and/or better control of the mem- 
brane's state. 
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Appendix A: Derivation of interaction Hamiltonian 

Our starting assumption is that the magnetic domain 
on the membrane is a point dipole located at the origin. 
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c. 3rd meas. F v = — 1, is = n/ijj m d. 4th meas. = l,t 4 = 7r/w ; 



FIG. 4: (Color online) Post-measurement Wigner distribution 
functions of the membrane, initially in a coherent state a = 
1 + il. Wc notice that repeated measurements do indead lead 
to some interesting non-classical states. The measurement 
results and evolution time before measurement are given for 
each plot. 



The magnetic field due to that dipole is 

B ( r ) = t- ~5 (3(Mm-f)f - Mn 
47r r J 



(Al) 



It is the small inhomogeneity of the resulting magnetic 
field along the long axis z of the condensate, and the 
resulting variation in Larmor precession frequency, that 
permits us to characterize the center-of-mass mode of 
oscillation of the membrane. 

For a dipole polarized along the z-axis (/i m = /i m z), 
at distance r = (x 2 + y 2 + z 2 ) 1 / 2 , the components of the 



magnetic field (Al) are 



B x (x,y,z) = 


MO Mm 

4tt r 3 


2>xz 


1 


(A2) 


B y (x,y,z) = 


Mo Mm 


'3yz~ 




(A3) 


4tt r 3 




1 


B z (x,y,z) = 


Mo Mm 
4tt r 3 


~'3z 2 


- 1 . 


(A4) 



Let x = Xq + x m , where Xq is the equilibrium value of 
x for a condensate atom, x = is the equilibrium po- 
sition of the membrane, and x m is the small sinusoidal 
displacement of the membrane around the origin. Ex- 
panding the expression for magnetic fields for small x m 
(up to first order) we get 



B,, 



B, 



Mo Mm 
47rr^ 

Mo Mm 

47T7q 

Mo Mm 
47rrfj 



3xqz 
[(2z 2 - 



3z(4xl 



Ibxoyz 



y 



z 2 ) 



(A5) 

(A6) 
(A7) 



3x (xl+y 2 -4z 2 ) 



Here, tq = (xQ+y 2 +Z 2 ) 1 / 2 . Assuming a two-dimensional 
condensate (i.e. y s» 0), we can set B y — > 0. For theo- 
retical simplicity, we assume that the measurements are 
made on the part of the condensate directly above the 
magnet, i.e. close to z = 0. Under these approximations, 
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the magnetic field simplifies to being predominantly in 
the z-direction, 

B z *^[~x + 3x m ]. (A8) 

Combining this with the quantization field B , we end 
up with the interaction Hamiltonian of the form 



V = ^b9fF z 



Bn- 



We can break this total magnetic field into a constant 
part, B c (sum of the first two terms), and a component 
that varies at the membrane frequency, B' v x m . Here, the 
gradient of the magnetic field is 



B' = 



3flQfI rn 



(A10) 



(A9) This is the form used in Eq. Q. 



Appendix B: Derivation of Successive Measurement Density Matrix 



Here we outline the steps used in section IV to derive Eq. (23). We arrive at the result via induction. Suppose that 
n — 1 measurements have been performed already. We denote the density matrix at this moment by p^"" 1 - 1 . After the 
next period of free evolution and the nth measurement, we have 

p (n) = M^U(t n )p( n -Vu\t n )MW. (Bl) 
To find the matrix elements, we insert completeness relationships 

(a,x f \pM\p,Xi) = J dx J dy(a,x f \M ( ^U(t n )\a n ,x)(a n ,x\p^\/3 n ,y)((3^ 



Now, the propagator can be derived from Eq. ( 12 ) by multiplying by the trace, substituting 77 



iu> m t/2, and 

including the contributions of the independent parts of the Hamiltonian. Using these techniques, we arrive at the 
needed matrix elements: 



(a,u\M^U(t n )\r,v) = M^/ . f?"" exp 

tirhsmLu m t 



t n n L0 T - t n 5VLT 2 



+ 



Ah 



(u + v + 2At) 2 tan ^ - (u - v) 2 cot ^ 



(B3) 



The fact that Eq. (12) and the individual terms of Eq. (14) (i.e. terms of the form of Eq. (13)) do not differ except 



for shifts in the coordinates and some complex phase shifts leads us to conjecture that additional evolution and 
measurements will not alter the underlying structure of the density matrix. Thus, we try a matrix of the form ( 23 ) 
far//™ -1 ), but with unknowns for X^ 1 ^ 1 ^, etc. (superscript added for clarity) and see if evaluation of the integral 
(B2) produces a similar output for p( n '. Indeed it does, and it also gives us the recurrence relations used to derive 
Eqs.( [25p7l ). These are 



X in) [a] = X {n - 1] [a] cos uj m t n + P (n - 1] [a] sinw m <„ + ct„(1 - cos u) m t n ), 

pW[a] = P^-^ [a] cos u; m t n - A^-^sinw™^ + CT„sinw m i n , 

(/,<.") [a] = ^ n -^[a] + (X^l [a] 2 + X^-^ia} 2 ) cot u: m t n - 2A<"> [a]X (n - 1 '> [a] esc w m t n . 



(B4) 
(B5) 
(B6) 



The solutions to these recurrence relationships are the equations given above. 



Appendix C: Derivation of Coherent State Evolution 



We again will arrive at the result ( 34 ) by using induction and finding a set of recurrence relations as in Appendix 
B, but the steps involved are slightly different. We first split the Hamiltonian ^ into two parts Hjsec arL d (H m + V). 
These commute, so the unitary evolution operator can simply be factorized into two operators, one acting only on the 
BEC and another acting as a shifted harmonic oscillator. We also note that the unperturbed membrane Hamiltonian 
and H m + V are related by a simple unitary transformation, namely 

B - + * = D (-ir/') H - D {^/-)' (C1) 
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Using this fact and the well-known results 

exp[-%H m t/h]D(Q = Die-^Qexpl-iHnt/h], (C2) 
D(C)D(0 = D(C + Oe (a '- Ci)/2 , (C3) 



we can postulate that the system after n — 1 measurements is specified by a superposition of states as in ( 34 ) , apply 
an additional evolution and measurement, 

|*(»)) = M^U BEC (t n )U m (t n )\¥ n -V), (C4) 

and then use the result to find the following recurrences: 

a («)[a] = a^-^ia] cos io m t n + b^-^[a] sin ij m t n - ^ (1 - cos w m t„), (C5) 
b (n) [a] = 6 ( " _1) [cT]cosa; m f n -a (n_1) [o-]sinw m t„ - ^— -wau) m t n , (C6) 

0(n) = Q (n-1) +Qn ( 6 (n-1) ( C7 ) 

Note the strong resemblance between the first two of these equations and those from the previous appendix. This is 
to be expected, as the real and imaginary parts of the displacement are proportional to the expected position and 
momentum of the coherent state, respectively. The solutions to these equations are Eqs. ( 35p7 ). 



Appendix D: Experimental Implementation 

Here, we provide details on the experimental detection 
of the quantum backaction on the micromechanical mem- 
brane. The membrane is composed of Silicon Nitride with 
a micron-scale magnetic domain deposited at its center. 
These structures are fabricated by coating silicon wafers 
with around 100 nm of high stress silicon nitride and then 
etching the silicon in a 100 /mi x 100 /im region to reveal 
a free standing membrane. Quality factors exceeding 10 6 
have been demonstrated with similar SiN membranes for 
resonance frequencies in the range of 500 kHz [18 . 

The magnetic domain is deposited on this membrane 
by a combination of photolithography and sputtering of a 
high permeability material like permalloy |17| . Peak sur- 
face fields on the order of 1 T at the surface of this mag- 
netic domain should ensure field gradients around 10 3 
T/m in the regions above the domain. The membrane 
is supported on a cryogenically cooled flange housed in 
a UHV chamber. Sputtered gold films deposited on the 
substrate in the regions surrounding the membrane serve 
as mirrors for trapping and imaging. Spinor Bose con- 
densates of 87 Rb in the F — 1 hyperfine state will be 
optically confined about 5 /im above the membrane in 
a quasi-2D trap created by shallow angle interference 
of laser beams reflected off the gold surface [TI5]. This 
ensures a quasi-2D confinement of the condensates such 
that the spatial extent of the gas in the direction nor- 
mal to the membrane is less than the spin healing length 
£ = yJ%l1m\c<2Xi\. Here, C2 = 47rfi 2 (<i2 — ao)/m is the spin 
dependent coupling strength of the spinor gas and 00(02) 
is the s-wave scattering length in the F — 0(F — 2) chan- 
nel. In this way, we can freeze out spin dynamics along 



this dimension [12"] . 

The estimation of the membrane's position is based on 
detecting the Zeeman shift at the location of the conden- 
sate due to the micromotion of the membrane. This is 
achieved by magnetization-sensitive phase contrast imag- 
ing as demonstrated in |12j . To summarize this technique 
briefly, a sequence of non-destructive phase contrast im- 
ages of the condensate are obtained with far off-resonant 
circular polarized light. Due to Larmor precession of the 
gas, the phase imprinted on the probe light is modulated 
at the Larmor frequency leading to a phase contrast sig- 
nal given by 

s = l + 2na ( 7 /2A)[a + ai (F v ) + a 2 {F*)} (Dl) 

where n is the column density of the gas, Co = 3A 2 /2-7T is 
the resonant cross section, A is the detuning of the probe 
light from resonance and 7 is the natural linewidth. F y is 
the projection of the local atomic spin along the imaging 
axis. The constants ag, a± and 02 depend on the detuning 
of the probe light and describe the isotropic polarization 
and optical activity of the condensate. By analyzing the 
sequence of phase contrast images on a pixel-by-pixel ba- 
sis, we can estimate the local Larmor precession rate and 
hence, the local magnetic field. At the shot noise limit, 
we estimate a magnetic field sensitivity on the order of 
15 fT/Hz 1 / 2 for a condensate of area 50 /mi 2 in the plane 
of the membrane [20]. For the expected field gradient of 
10 3 T/m at the location of the condensate, this field sen- 
sitivity translates to a position sensitivity of 15 x 10~ 18 
m/Hz 1 / 2 . This is much smaller than the zero point mo- 
tion of the oscillator which has an amplitude on the order 
of 6 x 10" 15 m. 
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